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Abstract 

The damping rate and free path of neutrinos in the singlet Majoron model 
have been calculated including both finite temperature and symmetry breaking 
effects. The behaviour of right- and left-handed fermions are found inherently 
different. While the damping rates of the left-handed leptons are essentially 
model independent, e.g. directly applicable to the Standard Model, for the 
right-handed particles the rates are crucially sensitive to parameters of the scalar 
sector. In general, the damping rates are fairly large. The possibility of the right- 
handed neutrinos to penetrate deep into the broken phase in the electroweak 
phase transition still remains, however, for some parts of the parameter space. 
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1 Introduction 



A particle traversing a thermal medium undergoes continual incoherent scatterings 
with plasma particles. Because of these interactions energy and momenta of the par- 
ticle are not sharply defined but have a finite "width", 7. The quantity 7, called 
the damping rate, describes the "decay" of the state (quasiparticle) into the plasma. 
Only when 7 is small compared with the energy will the quasiparticle propagate over 
meaningful distances in the plasma as a true physical excitation and with significant 
physical effects. 

The calculation of the damping rate of quasiparticle at zero spatial momentum, 
7(5cd(0), in QCD plasma has been a widely disputed issue in literature [|l], @]. The 
controversies in the results obtained for this quantity were solved by Pisarski and 
Braaten and Pisarski P, ^ who developed a general method for resumming the higher 
loop diagrams which contribute to the damping rate in leading order. 

The method of Braaten and Pisarski can be applied not only to QCD but also to 
other interactions. In this article we shall use it to calculate the damping rates of low 
momenta neutrinos in the electroweak plasma of the early universe. 

We are interested in the determination of the damping rates in connection to 
the so-called charge transport mechanism for generating the baryon asymmetry of 
the universe during the electroweak phase transition 0. This scenario is based on 
the assumption that the electroweak phase transition is of first order. Bubbles of 
the broken phase are created by fluctuations, and they grow and fill the pre-existing 
unbroken phase. As a bubble expands a flux of leptons from the unbroken phase hits 
the wall separating the phases, and some portion of the flux is reflected. As a result of 
violation of lepton number and C- and CP-parities, a net number of leptons is created 
in the unbroken phase. This is turned to a net baryon number by the sphaleron 
transitions. The produced baryons pass to the broken region, where baryon number is 
conserved, resulting in non-equal numbers of baryons and antibaryons in the universe 
we live. 

The reflection is generally most efficient at small spatial momenta (that is why 
the damping rate at zero momentum is a relevant quantity to study) for the following 
reason. When a quasiparticle penetrates to the broken phase region it has an energy 
determined by the symmetric phase dispersion relation, until it begins to see the 
thermal background of the broken phase. If a particle has energy which is smaller 
than the thermal mass it acquires in the broken phase, it can propagate there only 
the time determined by the thermalization rate, after which it is reflected efficiently. 
Thermal masses are generally of the order gT, i.e. particles moving with small spatial 
momenta are reflected efficiently. 

If its coherence length ~ I/27 is small compared with its effective reflection length. 



a typical quasiparticle decays into the plasma before it is reflected. This would suppress 
the generation of the baryon number making the scenario less viable. Actually, it has 
been claimed |^ that in the picture where the baryon number is created, instead 
of neutrinos, by a reflection of quarks from the phase boundary, this is what happens: 
decoherence effects crucially diminish the baryon number creation. It is interesting to 
investigate what is the situation with neutrinos. In the present paper we will calculate 
the damping rates for left-handed and right-handed neutrinos in the Majoron model. 
The consequences of the result with respect to the baryon number generation will be 
analyzed in a separate publication p. 

The organization of the paper is the following. In Section 2 we will introduce the 
model we are working with, the singlet Majoron model. In Section 3 we derive the 
effective vertices and propagators, and in Section 4 we apply them to the calculation 
of the zero momenta damping rates. The numerical results are presented in Section 5, 
and the conclusions drawn in Section 6. 



2 The singlet Majoron model 

For evaluating the neutrino damping rate one has to specify the model. We will use 
the singlet Majoron model 0, which is an extension of Standard Model containing 
an extra gauge singlet scalar field S, as well as gauge singlet right-handed neutrinos. 
The lagrangian of the model is the most general one with the requirement that it 
has a global B — L symmetry. Spontaneous breaking of this global symmetry by a 
non- vanishing vacuum expectation value (vev) of the singlet scalar, (S) = /, generates 
Majorana masses to the neutrinos. The see-saw mechanism can be implemented yield- 
ing three neutrinos which are very light compared with the other fermions of the model 
and three neutrinos which are very heavy. It is usually assumed, to avoid an ad hoc 
hierarchy among the energy scales, that / is comparable in size with the vev ($) = / 
of the ordinary doublet Higgs $ responsible on the electroweak symmetry breaking. 



The model can fulfill the Sakharov conditions ||10| for the creation of a cosmological 
asymmetry between baryons and antibaryons, as was demonstrated in 

The singlet Majoron model contains the following new Yukawa interactions not 
present in the Standard Model: 

Cy = IlXd^Nu + NrXmSNI (1) 

Here $ is the conjugate of the Higgs doublet, is a vector of left-handed lepton 
doublets and is a vector of right-handed singlet neutrinos. The Yukawa coupling 
matrices and Am are in general non-diagonal in flavor space. 

The first term of (|1]) gives rise to Dirac masses and the second term to Majorana 
masses for neutrinos. In order to satisfy the experimental mass limits of the known 



neutrinos one has to require, given our assumption that there is no strong hierarchy 
between the vevs / and /, that the elements oi Xd are much smaller than the elements 
of Am, so as to make the predominantly left-handed neutrinos very light. A detailed 
discussion of the masses, mixings and phase transitions in the singlet Majoron model 



can be found in [TTJ. We just note that due to a large number of free parameters 
appearing in the model and a lack of constraining experimental data, the vacuum 
expectation value /, singlet scalar mass and Majorana masses can be considered as 
practically free parameters. 

3 Effective vertices and propagators 

The damping rate is defined as the imaginary part of the solution u = uj{k) of the 
dispersion relation 

det {^j\uj, k)) = det (cJ7° - k ■ 7 - k) - m/) = 0, (2) 

where Aj is the fermion propagator, S is the fermion self-energy and mj is the vacuum 
mass of the fermion. The contributions to the finite temperature self-energies that 
one has to consider in order to determine the zero momentum damping rate of the 
Majorana neutrinos in the singlet Majoron model are displayed in Fig. 1. In this 
section we shall outline the computation of the effective propagators and vertices 
(presented by the 'blobs') out of which these amplitudes are constructed. In the next 
section we will apply the results to evaluate the damping rates. In our calculation we 
will make the following approximations. We will compute all quantities in leading order 
in coupling constants, neglect the Yukawa couplings except the Majorana couplings of 
the right-handed neutrinos, and when necessary take subleading terms into account 
to avoid problems with infrared divergent graphs. 

We will apply the method of Braaten and Pisarski |p to identify those higher- 
loop diagrams, called hard thermal loops (HTL), which contribute in leading order 
and which are necessary for the gauge invariance of the damping rate. By definition, 
HTL's are those diagrams which contribution to effective vertices is of the same order 
in coupling constants as the contribution of the tree level diagrams. According to @], 
an amplitude contributing the effective vertices and propagators can be an HTL only 
when its external legs are all soft, i.e. the momenta of the external particles are of the 
order ~ gT, where T is temperature and g is a generic coupling constant. 

Let P = {po, p) denote the external soft momentum and K = (fco, k) the loop 
momentum. Whether a given amplitude is an HTL can be determined by applying the 
power counting rules introduced in Ref. [^: the integration element / d^k contributes 
T^; the first propagator in the loop and the summation over ko together contribute 
1/T, and all other propagators give a contribution 1/(PT) each; in the numerator 



contributes T; in the numerator contributes gT; for loop where internal lines are 
all bosonic or fermionic there is an extra factor of P/T. 

The hard part of an amplitude corresponds to making a linear approximation in 
the external momentum P ~ gT in propagators, i.e. neglecting P^ as small compared 
with P ■ K. The particle masses, which appear in the propagators quadratically, can 
be neglected, because they are generally of the order gv, where f is a generic vacuum 
expectation value of a symmetry breaking Higgs boson and g is a, generic coupling. 
An exception seems to be scalar masses, which are lower order in coupling constant, 
m| ~ gv"^. However, perturbatively in finite temperature one has v ~ g^^'^T, and 
since the soft scale of scalar lines should be identified with ~ g^^'^T, we can neglect 
in leading order also the scalar masses in loops. For example, calculating the scalar 
tadpole diagram, which arises from scalar self-coupling, gives for the scalar a thermal 
mass ~ g^^'^T, which is small compared with the relevant loop momentum. 

Let us now move to consider in one loop level the effective two-point, three-point 
and four-point functions entering the diagrams of Fig. 1. There is obviously a great 
number of loop diagrams contributing, but it turns out that just a handful of them 
are HTL's and thus relevant in leading order evaluation of the damping rates. We 
will not go through our analysis in full detail diagram by diagram here, but rather 
just describe the general methods we have used to discard the amplitudes yielding a 
non-leading contribution and to pick up the HTL's. 

Let us first consider diagrams for the left-handed neutrinos. In the case of the gauge 
boson - neutrino pair three-point function there are five different types of diagrams, 
which may be conveniently classified according to the types of their internal lines. Let 
F denote a fermion line, S a scalar line and A a gauge boson line. The diagrams 
with configurations FFA and FAA have good infrared behaviour. They are HTL's, 
and the calculation of them is presented in the Appendix. The diagrams FFS and 
FSS have also good infrared behaviour and they are HTL's, but because they contain 
small Yukawa couplings we neglect them. The configuration FAS has a vertex V(AAS) 
which originates in the symmetry breaking and which has with a coupling proportional 
to g^f. By the power counting rules its contribution would be ~ g'^f/T, i.e. of higher 
order than the tree level contribution. However, due to the behaviour of this graph in 
the infrared region, its contribution is larger than given by the pure power counting 
argument. To see this, and to introduce our scheme to approximate the contributions 
from the infrared integration regions, let us consider the amplitude FAS in more detail. 



Using the standard imaginary time formalism |I2| to perform the summation over 
discrete loop energy (see also Appendix), the diagram FAS may be cast into the form 



r(FAS) ~ g'f I 



(2^ 



' I fBjuJk) 1 1 

(3) 



Here = (k + p)^ + m^, fs is the Bose- Einstein distribution function, the Pj's are 
some combinations of the external momenta, and m is the mass of the loop particle. 
The external energies p° are here continued to real values. If one makes a linear 
approximation in the external four-momenta and neglects the masses in denominators, 
the integral in (|^) becomes infrared divergent and it behaves as 

r,FAS),„^///f|;^4/f, (4) 

where P denotes generically the soft external momentum scale. As the integral is only 
logarithmically divergent the contribution from the infrared region is of the order of 
\og{g)f/T. Hence r{FAS) is not an HTL in spite of its divergent infrared behaviour. 

By generalizing previous example one obtains power counting rules for the in- 
frared region: the summation over ko and one propagator contributes together treating 
l/{gT)] the integration element d'^k contributes (gT)^; the loop momentum K is taken 
to be of the order of the external soft momentum P ~ gT both in the nominator and 
in the denominator; and the Bose-Einstein distribution function contributes 1/g. 

There are five types of diagrams contributing to the scalar-fermion pair three-point 
function: FFA, FFS, FAS, FAA and FSS. All but one, FAA, of these configurations 
contain small Yukawa couplings, and hence we shall neglect them. The diagram FAA 
is due to the symmetry breaking, and it behaves as ~ g'^f/T both in the ultraviolet 
and infrared regions. Because also the tree level coupling of the left-handed neutrinos 
and a scalar is proportional to small Yukawa coupling, the effective scalar - neutrino 
pair vertex can be omitted for the left-handed neutrinos altogether. 

A note concerning the configuration FSS is in order. It contains a scalar three 
point vertex arising through the symmetry breaking from the scalar four-point cou- 
pling, generically denoted by A, of the symmetric theory. Power counting results in a 
contribution ~ gv/T , where w is a generic vacuum expectation value parameter. Thus 
if V is of the same order than temperature T, one should consider the FSS diagram as 
a HTL. But here one should be careful to check whether the contribution comes from 
the scalar or from the fermion line (see Eq. ( |A10| )). By identifying the soft scale of the 
scalars with A^^^T, these contributions are ~ (/yf/T and ~ X^/'^gyv/T^ respectively. 
Here gy is a generic Yukawa coupling. One can check that also the ultraviolet region 
gives a subleading contribution and hence FSS diagram is not an HTL. Formally the 
same result is achieved if A^^^ is considered as the perturbation parameter of the scalar 
sector. 

With the above note in mind one can check that from the 13 possible diagrams 
of the two gauge boson - neutrino pair four-point function only three, FAAA, FFAA 
and FFFA, contribute. The calculation of these is outlined in the Appendix. There 
are five diagrams which have vertices arising from the symmetry breaking: FFAS; 
two diagrams of the type FASS; two diagrams of the type FAAS (these are neglected 



because they are not HTL's); two non-HTL diagrams with vertices of the symmetric 
theory: FAA and FSS; and three diagrams containing small Yukawa couplings: FSSS, 
FFSS and FFFS. 

Further one can check that none of the 13 possible diagrams of the two scalar - 
neutrino pair diagrams contribute. There are six diagrams with symmetry breaking 
vertices: FSSS, FAAS, FAAA, FASS, FFSS and FFAA; these are again neglected 
because their are not HTL's; three non-HTL diagrams of the symmetric theory: FAAS, 
FSS and FAA; and four diagrams of the symmetric theory containing small Yukawa 
couplings: FASS, FFAS, FFFS and FFFA. 

Finally, none of the three possible diagrams FFSS, FFAS and FFAA contributes to 
the effective four neutrino vertex . Thus we are left with the effective contributions 
to the self-energy of the left-handed neutrino of Fig. 1 a. 

Similar considerations of the effective vertices entering to the calculation of the 
self-energy of the right-handed neutrino are much more simpler because the particle 
is a gauge singlet: there are no three-point or four-point functions of right-handed 
neutrino pair and gauge bosons. 

The three-point function of a right-handed neutrino pair and a scalar has two 
possible one loop contributions: the diagrams FFS and FSS. The latter is not HTL 
by the argument applied already in the left-handed case. The former is neither HTL: 
the fermion propagators produce a term to the denominator of the integrand of 
the diagram and hence cancel corresponding factor in the nominator. Effectively this 
gives an extra factor of g to the order this diagram. Thus the effective three point 
function of the right-handed neutrino pair and a scalar consists just of the tree level 
term. 

The four-point function of a right-handed neutrino pair and two scalars has four 
possible one loop contributions: the diagrams FSSS, FSS, FFSS and FFFS. The first 
two of these are non-HTL diagrams by power counting, while the third and the fourth 
are not HTL diagrams because of the term in the denominator due to the fermion 
propagators. 

Finally, similarly as in the case of the case of left-handed neutrinos the possible 
contribution FFSS to the four neutrino vertex is not a HTL. Thus we are left with 
the effective contribution to the self-energy of the right-handed neutrino of Fig. 1 b. 
The effective neutrino pair - scalar three point function consists only of the tree level 
contribution. 

The two-point functions which contribute to the effective propagators of left- and 
right-handed neutrinos, scalars and gauge bosons are selected correspondingly. It turns 
out that only the diagrams of the symmetric theory contribute in leading order. Fur- 
ther, in the case of left-handed fermions, because of the negligible Yukawa couplings, 
only the diagram of the type FA contributes and, in the case of right-handed neu- 



trinos, being gauge singlets, only the diagram FS contributes. In the case of scalars 
the two-point function is independent of the external momentum and it just gives a 
contribution to the mass term of the scalars. 

The effective, high-temperature, two-point functions of the gauge bosons and fermions 
are well known (see e.g. Refs. [0, |l^). They are listed in the Appendix. 

4 Zero momentum damping rate 

The neutrino damping rates at zero momentum are calculated to leading order by 
evaluating the one-loop diagrams of Figs. 1 a and 1 b, where the effective couplings 
and propagators include the hard thermal loops discussed in the previous Section. 
We will consider the damping rates of the left-handed neutrinos and right-handed 
neutrinos, that is, of the light and heavy Majorana neutrinos, separately. There is an 
essential difference between these two cases because the interactions of the two sets of 
particles differ. 

4.1 Left-handed neutrinos 

In the left-handed sector, according to the analysis of the previous Section, only the 
graphs of unbroken gauge theory contribute. Thus the calculation of effective n-point 
functions proceeds similarly as in the case of QCD except that the group theoretical 
factors are replaced by those corresponding to the SU{2) xU{l) symmetry, and one has 
to take into account the left-handed nature of neutrinos, as well as the scalar doublet 
contribution to the gauge boson thermal masses. Although the modifications to the 
QCD results are quite straightforward, we present for completeness some technical 
details in Appendix. 

We wish to study the leading symmetry breaking effects, that is, an inclusion of 
the vacuum mass terms of gauge bosons, to the zero momentum damping rate of light 
neutrinos. Although the diagrams arising from the spontaneous symmetry breaking 
were neglected as subleading when the effective vertices were calculated, the gauge 
boson and fermion masses created through the symmetry breaking, so-called vacuum 
masses, may have as important effects to the damping rate as the thermal masses 
arising from the interactions with the plasma. One should note that the inclusion of the 
vacuum mass terms does not spoil the gauge invariance of the damping rates, since the 
Ward identities obeyed by the n-point functions remain unaltered [H . These identities 
guarantee that the gauge dependent piece of the gauge boson propagator gives zero 
contribution to the on-shell damping rate. The relevant identities are displayed in 
Appendix. 



The effective self-energy of a left-handed neutrino may be written as 



s(i)(p) = / tSf P - K)MP - K)rf (P, P - K)GZiK) 



+ \^t^{P,K)G^j:B ■ (5) 



Here is the effective gauge boson - neutrino pair three-point function, F^^f is the 
effective four-point vertex for neutrinos and gauge bosons, A; is the effective propa- 
gator of a left-handed neutrino, and G^^^ is the effective propagator of gauge bosons. 
The indices A, B refer to the SU{2) x U{1) gauge bosons. 

The gauge boson propagator may be divided into longitudinal and transverse pro- 
jections ||r3| : 

G^^^{K)=^^^P,^ + /\i^Q,^, (6) 

where P^^, Q^i, together with K^^ = k^ky/K"^ form a set of projection operators 
obeying P^y + Q^^ + K^^ = g^^. Explicitly, 

Pf.. = 9,.. - U^Uy + ^{K^ - koU^){Ky - koUy), (7) 

where U = (1,0,0,0) is the four-velocity of the plasma. The gauge dependent part, 
which is proportional to K^y, is not displayed in Eq. (|]). The longitudinal and trans- 
verse propagators Al,t have the form 

At%{k) = {K'-M'-U,AK)y^l, (8) 

where is the vacuum mass matrix, and IIl,t are the longitudinal and transverse 
projections of the vacuum polarization tensor. 

The damping rate is defined as the imaginary part of the effective self-energy, 
which one finds as the solution of the dispersion relation. It reads at zero momentum 
as follows: 

7(0) = ^Tr(7olmS«(po = ^«,P = 0)) 

= ^DiscTr(7oS«(po = cc^^,P = 0)). (9) 

The operation Disc/(x) = f{x + iO^) — f{x + iO~) extracts the disconnected part of 
the function /, and ui is the thermal mass of the left-handed neutrino. After some 
tedious algebra, we may write the damping rate to the form 



1 „^ r d^k 



7(0) = - Disc / 



(2vr) 



(10) 



where Fa and Fb are group generators (multiplied with the corresponding gauge cou- 
pling constant) and 



^ Y^^^{P-K)^^\{po-k,f-{k±uj{) 



A;2 



In 



The neutrino propagator is here divided into the helicity components Aj., 
AKP) = ^A+(7° - p • 7) + ^A_(7° + P • 7). 



'Ill 



(12) 



The form Eq. (|T0|) is similar to the QCD case 0, only that the gauge boson propagator 
factors have now a more complicated form. 

The sum over the discrete loop energy values may be performed with the help of 
the formula 

DiscT^i?(A:o)F(po-A:o)|,„.K = 

fco 

dujfB{uj)pB{uj) / duj'fF{uj')pF{uj')6{E - iO - Uj'), (13) 

-OO J — oo 

where the argument of the function B (F) is even (odd) multiple of uTi, fB{F) is the 
Bose-Einstein (Fermi-Dirac) distribution, Pb{^^) = l/{2'n'i)Disc B{uj) is the spectral 
function of B{uj) and pp is similarly the spectral function of F{u!). 

The damping rate of the left-handed neutrino in the zero momentum limit is thus 
given by the formula 



iovr 



(14) 



where 



LUt 



UJ 



k dkdLjdLj'6{uJi — UJ — uj')fp{uj')—fB{uj) x 



2pl{uj, k) 



UJ 



J2p±{^',k) {uji+uj'Tkf + 



Pt{uj, k) 



j:p±iuj',k)^{uj"-{k±uj,ry+ 



UJ 

2 e - 



and 
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PT 
P± 



fc2 

Disc 
Disc 
Disc 



^^{FA.FB}^r{uj,k) 



k^ 

UJ^-k^ 



ufeik'-uj") 



(15) 



k^ 



^^{FA,FB}A^^{uj,k) 
[A±iuj',k)]. 



(16) 



Here g2 is the SU{2) gauge coupling.Note that Eq. (p^ apphes for the both com- 
ponents of the left-handed lepton doublet, the neutrino and the charged lepton, by 
choosing the appropriate T3 quantum number. 

The gauge boson spectral functions separate naturally to the W-boson, Z-boson 
and photon terms: 



(17) 



^ K='y,Z,W 



where the delta-function terms arise from the poles of the propagators Al^t and the 
step function term from the cut fc^ > tu^ of the logarithmic terms in the propagators. 
The explicit expressions of these are straightforward to calculate, resulting, however, 
in rather lengthy formulas not displayed here. Similarly, the fermion spectral functions 
p± may be expressed as 

p±{uj, k) = Z+{k)6{uj T io+{k)) + Z.6{uj ± co^k) + e{k'^ - uj^)R±{uj, k) (18) 

with R_{ijj, k) = R+{—Lj, k). 

Following Ref. |]r3|, the residues ^2'^^^ and Z± may be interpreted as proba- 
bilities of creating an excitation obeying the energy-momentum dispersion relation 
u = uj]^^'^ (k) and u = ±uj±{k) in the plasma, respectively. In the case of fermions 
we may interpret the two branches ±uj± as particle-like and hole-like excitations. The 
negative-energy solutions are interpreted as antiparticles as usual. Further, the coeffi- 
cient of the step function may be interpreted as the contribution of the particles in the 
heat bath to the spectral functions: the quasiparticles are not strictly on mass shell, 
because they interact with the heat bath. 

With these interpretations, one may divide the coefficient a into separate contribu- 
tions according to whether the energies u and uj' in Eq. (p!5| ) are off-shell or on-shell. 
Particularly, if the both energies are on-shell, the corresponding contributions are in- 
terpreted as processes Uplh ^ A and Up Ip — > A, where A denotes generically one of the 
bosons 7, Z or W , I denotes generically the left-handed doublet lepton, and subscripts 
p and h refer to particle and hole excitations, respectively [|T6|. The other possible 
processes are kinematically forbidden, and it is easy to show that the processes above 
are allowed only if the thermal mass of the boson is larger than 2uji. Note that the hole 
excitation has same quantum numbers as the anti-particle excitation but the opposite 
helicity [|T5[] . 



In addition to the on-shell contribution there are contributions where at least one of 
the particles in the loop in Fig. 1 a is off-shell. These contributions may be interpreted 
as due to the interactions of the quasiparticles with the heat bath. The behaviour of 
the integrand in Eq. (|I5]) for these contributions is such that the main contribution 
to a comes from the soft region. In that case one can approximate exp(co'//T) ~ 1, 



fsiuj) ~ T/cu and = 1/2- The result is that the parameter a is indepen- 

dent of temperature (the same is the case in QCD The on-shell contribution is 
exceptionah schematically it reads 

«os = f^'Tff^^l^i /H-^^(fco))/BK,T(fco)) X (...), (19) 

where fco is the solution to the equation uji + uj±{k) —oji^Tik) = 0, prime denotes deriva- 
tion and the ellipsis represent contributions that are not relevant in our reasoning. It 
can be shown that and approach unity when k increases, while Z^ and Z_ 
approach zero exponentially. Further, ujt and 00+ approach the free particle dispersion 
relation lu = (fc^ -|- M^)^/^, M is the vacuum mass, with increasing k, while lul and 
LV- approach the massless dispersion relation uj = k exponentially. Consider the case 
of transverse boson and particle-like neutrino excitations in Eq. ([T9|). If the vacuum 
mass of the gauge boson is large, M ^ ui, then ko ^ M"^ /{2uJi) and 



2 



aos ~ ^e-'^"/^, (20) 



which means that without Boltzmann suppression a would grow to high values with 
increasing / ~ M. 



4.2 Right-handed neutrinos 

The calculation of the zero momentum damping rates of right-handed neutrinos pro- 
ceeds along the same lines as of the left-handed ones, the main difference being the 
inclusion of the vacuum mass term to the propagator of the right-handed neutrino. 
According to the analysis of the previous section we need to calculate the disconnected 
part of the diagram in Fig. 1 b. In the case of right-handed neutrinos the three-point 
vertex is bare; there is no HTL-contribution to it. Furthermore, as mentioned in the 
previous Section, the HTL-contributions to the scalar self-energies are independent of 
the external soft momentum P. They merely add a contribution to the mass term 
of scalars, which is, in fact, already taken into account when we use the temperature 
corrected scalar potential. 

The evaluation of the effective self-energy of the right-handed neutrinos follows 
the same lines as that of the left-handed neutrinos. The only change is the replacement 
of the squared thermal masses with 

= AmAL^. (21) 



Because we assume that the Dirac mass terms connecting the left- and right-handed 
neutrinos are neglected as compared with the Majorana mass terms, the vacuum mass 



terms and the thermal energies are diagonahzed simultaneously. The effective la- 
grangian can then be presented in the mass eigenstate basis as follows: 

^e// = \xi{l-P- Sl) XI + \xh{l-P-^R- mu) Xh + --- (22) 

Here xi Xh represent the light (mainly left-handed) and heavy (mainly right- 
handed) mass eigenstates of neutrinos, respectively. The self-energies and are 
diagonal. Thermal masses of the the heavy flavour eigenstate neutrinos are (i=l,2,3) 

= (23) 

Here mM,i is the vacuum mass given by -\/2/A^^^|, where A^!f^ is the diagonalized form 
of the coupling matrix of the right-handed neutrinos and the singlet scalar defined in 
Eq. (0). 

The vacuum mass terms of the scalars give rise to a mixing between the singlet 
Higgs and the neutral member of the Higgs doublet. Hence the scalar propagator 
appearing in the neutrino self-energy amplitudes is a 2 x 2 matrix. The component of 
this matrix corresponding to the dominant channel to which heavy neutrinos couple 
with a large Yukawa coupling is given by 

^24) 



Here V^h = —m\f^{T) and = K"^ — m1g{T) are the inverse propagators of scalars 
h and s, and T^s is the mass mixing term arising from the scalar potential through 
the symmetry breaking. Subscripts h and s refer to the neutral, P-even components 
of the doublet and singlet fields, respectively. As a matter of fact, one can neglect in 
( p4D the scalar mass mixing term F/j^. As will be seen in the next Section, inclusion of 
it would play subdominant role as compared with the much more substantial effects of 
the inclusion of the vacuum mass term of the right-handed neutrino. In addition to the 
scalar contribution in Eq. (p^), there is also a Majoron contribution. We shall neglect 
the thermal mass differences of the Majoron and the other scalars, and combine the 
effects of these particles by using a complex scalar field S to describe the scalar degrees 
of freedom. 

Unlike for the left-handed neutrinos, the zero momentum damping rates for the 
particle- and hole-like excitations differ in the case of the right-handed neutrinos as a 
result of the inclusion of the vacuum mass term. Namely, the imaginary part of the 
self-energy at zero spatial momentum has two contributions ai and ^2 defined 
through 

ImS^')(po = ujf, p = 0) = + 4, (25) 



where 



(26) 



are the thermal masses of the particle (+ sign) and the hole (- sign) excitations in the 
case of nonvanishing vacuum mass. From the zero momentum dispersion relation 

Po lai] - [mM + ^012} = (27) 

one may then deduce that the zero momentum damping rates 7±(0) of the right-handed 
neutrinos read 



7±(0) 



("l ±"2 



f28) 



Similarly as in the left-handed calculation the damping rates can be written in the 
form 



7±(0) 



lA 



diag 1 2 
M 



T 



IGtt 



-a±, 



(29) 



where 



4(e^'- + 1) X / edk / dujduj'5{LU^ -LU-Lu') 



I 00 

Here the scalar spectral function ps has a very simple form: 
ps{uj,k) 1 



2ul{k) 



{6{uo - ujs{k)) + 6{uo + uos{k))) 



(30) 



(31) 



where uJs{k) = {k'^ + ms{TYY^'^. The neutrino spectral function contains two parts, 
Pri and pr2, corresponding to the imaginary parts ai,2 in Eq. (^51) . The right-handed 
neutrino propagator may be written in the form 



(/(Po,p)7o + h{Po,p)P -l-ruM 
1 



(/(Po,p)7o + /i(Po,p)p-7 + "^M), (32) 



P{Po,p) - h^{Po,p) - 

where functions / and h are defined through Eqs. ([A23|) and ( [A24| ) . One can then 
write the spectral functions Pri and pr2 as follows: 



Prl{uJ,k) 
Pr2(t^, k) 



Disc 
"27r7 

Disc 



f{u;,k) 



k) — h'^{u, k) — m 

mM_ 

P{uj, k) — h'^{uj, k) — m 



(33) 



Factorization to the helicity components is not possible in the case of the right- 
handed neutrino propagator because of the vacuum mass term. However, as before 



the propagator has four poles oj = ±uj^{k), which may be interpreted as dispersion 
relations of the particle-, antiparticle-, hole- and antihole-like excitations, but in con- 
trast with the vanishing vacuum mass limit there are now two residues corresponding 
to the two contributions to the spectral functions p^i and Pr2- Similarly there are also 
two off-shell contributions proportional to the step function ^(/c^ — u"^). 

The on-shell channels contributing the damping rates are listed in Table 1. Simi- 
larly as in the calculation concerning the left-handed neutrinos, the on-shell contribu- 
tion to the parameter a± may be written schematically as 

«±,os = I . f^T+'\. M±ujHko))fB{±uJs{ko)) X (. . .), (34) 

where fco is the solution of the equation cj^ ±uj^{k) ±ujs{k) = 0. The sign choices de- 
pend on the channel in question. Again, the residue terms Zi 2 for particle excitations 
only are non- negligible for large k^. Further, the denominator in Eq. (|3^) approaches 
zero for large ko, if there is a relative minus sign between the two derivatives uj'^{ko) 
and uj'g{kQ). These two conditions are fulfilled in the case of the particle damping rate 
for the channel Upi/p S and in the case of the hole damping rate for the channels 
PhS Up and Uhi^p S. These channels give the largest contributions to the damping 
rates, as is evident from Figures 3 and 4. 

5 Numerical results 

The zero temperature damping rates left- and right-handed neutrinos derived in the 
previous section have quite complicated analytic form. In this section we analyze their 
properties numerically. Instead of the damping rates themselves we will present in our 
plots the coefficients a and a± defined in Eqs. (|l^ and (^). 

In Fig. 2 we have plotted the coefficient a, which is proportional to the damping 
rate of the left-handed neutrinos, as a function of the vacuum expectation value / of 
the doublet Higgs. As can be seen from the figure, a is essentially independent on /. 
This result can be understood due to the fact that, although the vacuum mass does 
suppress the damping rate, there is now a new thermalization channel open, which 
has no counterpart e.g. in the pure QCD case @]. The annihilation of the left-handed 
neutrino with a neutrino from the heat bath namely gives a large contribution to the 
rate. In fact, its contribution increases with increasing /, in contrast with all the 
other contributions which decrease. For //T ~ 1.4 it corresponds to about a half of 
the whole rate. For very larger values of /, however, also this annihilation contribution 
starts to diminish due to Boltzmann suppression. 

The corresponding results for the right-handed neutrinos are presented in Figs. 3 
and 4. The damping rate and thus the parameter coefficients a± depend in this case not 



only on the Yukawa coupling but also on the symmetry breaking mass ms^vac ~ Am/' 
as well as the on thermal mass ms^th ~ T of the singlet Higgs S, all of which can be 
viewed as free parameters. We have presented a+ and a_ as functions of the scalar 
mass m| = m| + m| for some representative values of the Yukawa coupling qXm 
and ratio r = rriM /iTiji{T), where m^j is the thermal mass and tum the vacuum mass 
of the right-handed neutrino. It is noteworthy that the ratio r does not depend on 
the Yukawa coupling constant but only on the ratio f/T. It varies in the range from 
r = to r ~ 6, the upper limit corresponding to its value at / ~ T. 

The damping rates of holes and particles behave quite differently when the pa- 
rameters are varied, the case of holes showing more structure. The general trend for 
holes is that the larger r, the larger is the damping rate. This can be seen in Fig. 
3 a. However, for each value of r there exists an interval in the values of the ratio 
ms{T)/mR{T) where the damping rate is suppressed. E.g., for r = 2 the damping 
rate is only ~ 0.1 when ms{T) /uijiiT) lies between values 2 and 3, while outside this 
range it is some two orders of magnitude larger. This is due to the fact that only soft 
processes are allowed in this specific parameter region. 

Also we find that the damping rate of the holes is quite sensitive to the value 
of the Yukawa coupling \m (Fig- 3 b). The particle damping rate instead depends 
essentially only on the ratio r. When r increases the damping rate decreases until a 
threshold mass is reached, as depicted in Figs. 4 a and 4 b. At that threshold the 
annihilation channel opens which drastically increases the interaction rate. However, 
the value of Am does not affect the damping rate but only marginally. Also for large 
r the damping rates for the particle-like excitation are generally smaller than for the 
holes. The main reason is that incoming particles tend to scatter off the holes, and the 
number density of the holes in the heat bath is suppressed compared with the number 
density of particles. 



6 Conclusions 

In the present paper we have studied the thermalization of low momentum neutrinos in 
the heat bath in the framework of the singlet Majoron model. This is a relevant issue 
in connection to the charge transport mechanism for generating the baryon asymmetry 
during the electroweak phase transition |^ when the reflection of neutrinos from the 
phase wall and from the matter in the broken phase is considered. 

We found that the damping rate of the left-handed neutrinos depends only on the 
temperature corrected vev / of the doublet Higgs field, not on other parameters of 
the model. This dependence is not particularly strong, the coefficient a defined in 



Eq. (|T^ varying between a = 6 and a = 3.5 when / increases from to 2T. The 



result a ~ 6, corresponding to the unbroken gauge symmetry, is quite close to the 



value of the corresponding parameter obtained in QCD for quark damping, a ~ 5.7 
(for three flavours) We found the damping rate of the left-handed neutrinos being 
determined solely by diagrams involving gauge bosons. Therefore our result can, up to 
some corrections C(l), be directly generalized to the left-handed charged leptons, too. 
Furthermore, while the result was derived in the Majoron model, it is valid for the 
left-handed leptons of the Standard Model as well, since the new degrees of freedom 
the Majoron model introduces play no role in the leading order calculation. 

The damping rates the right-handed neutrinos, which do not couple to the gauge 
fields, are determined by the Yukawa couplings and properties of the singlet scalar 
Higgs. Hence the damping rate is more model dependent than in the case of left- 
handed neutrinos. Depending on the values of various parameters, such as the vacuum 
mass and the vacuum expectation value of the singlet Higgs, one can find large as well 
as small damping rates for the right-handed neutrinos. 

The present analysis has been performed at the electroweak phase transition tem- 
perature T ~ 100 GeV, whereas the QCD calculations of Ref. was performed at 
QCD phase transition temperature T = 200 MeV Q. On the other hand, in the region 
where the gauge symmetry is broken, the mass of the weak gauge bosons tends to 
decrease the damping rate. The net effect is that the rough estimate for the thermal- 
ization time, I/7 ~ 0.1(GeV)~^, used e.g. in |]18|, is about correct for the left-handed 
leptons, while the thermalization time of the right-handed neutrinos may differ from 
this crucially. 

For each particle species the damping rate should be smaller than the particle mass. 
In the opposite case it would be meaningless to speak about particles traversing the 
medium. In the view were the particle mass is thought to be built up gradually as a 
result of sequential forward scatterings, particles would in this case decay faster that 
they are formed. For the right-handed neutrinos this means that the mass given by 
( p3D has to be smaller than the damping rate (^), implying Aa//o± < 47r. This cuts 
out some large values of a± as unphysical. 

A general conclusion of the present paper is that the damping is fairly strong, not 
only in QCD, but also in the electroweak sector. Hence the reflection probability of 
the neutrinos may in general be substantially reduced which may crucially affect the 
charge transport mechanism of the baryon asymmetry generation of the early universe. 
However, there is a region in the parameter space where the damping rate is suppressed 
allowing the right-handed neutrinos to penetrate deep into the broken phase. 



"'^It is noteworthy that if QCD calculations of Ref. had been performed at the electroweak 
phase transition temperature, the effective number of degrees of freedom in the heat bath would be 
increased. 
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Appendix 

In this Appendix we present the derivation of the effective vertices needed for the 
evaluation of the damping rate of the left-handed neutrinos. 

The hard thermal loops which enter into the calculation of the damping rate of the 
left-handed neutrinos arise from the sums of the form |17| 



Hfi'-''"-' =TtK^' ■■■K''"-'Y[i^{K + PiY -m'^ij , (Al) 

i=l 

where 

Tr = rV/-^. (A2) 

It is understood that the zero component of the momentum attached to a boson 
(fermion) line is an even (odd) multiple of ttTz, and the momenta Pi are linear com- 
binations of external momenta. 

Consider then the effective three-point function of an SU{2) gauge boson and a 
fermion pair. There are three loops contributing: one with two SU{2) gauge boson 
propagators and one fermion propagator; one with two fermion propagators and one 
SU(2) gauge boson propagator; and one with two fermion propagators and one U{1) 
gauge boson propagator. Extracting the HTL contribution and taking into account 
the group theoretic factors we obtain in the case of SU(2) gauge bosons 

r^'^ = g2T^ (g'" + 4C^ff3l2) lu, (A3) 



where is an SU(2) generator and Cp = 3/4(7| + Y'^/Agf is the Casimir operator 
of the (2, y) representation of SU{2) x U{1) scaled with the appropriate couplings. 
In H^'^ the first subscript indicates the number of internal lines and the second one 
denotes the number of fermion propagators. The property i?^.'^ = —Hl^'^n-m "W'as used 
in derivation of the Eq. (|A3|) . A similar formula is obtained for the effective three-point 
function Fi of a fermion pair and a U{1) gauge boson: 

T^,=gi\{g^'' + ^C'pmi)l.. (A4) 

Note that the three-point functions above are of the general form 

f:^ = Fa {g^" + 4C^i73^^) 7., (A5) 



where Fa denotes the generator of the group in question multiphed by the correspond- 
ing couphng. 

Consider then four-point function of two gauge bosons and a fermion pair. There 
are three diagrams (+ crossed channels), contributing: those with one, two and three 
fermion propagators. Extracting the HTL-contribution, calculating the group theo- 
retic factors and treating the SU(2) and f/(l) gauge bosons equally the four-point 
function may be cast to the form 

^AB = -8 {-[Fa, Fc][Fb, Fc] + Fc{Fa, Fb}Fc) i/^^f 7a 

+ 8 {[Fbi Fc\FaFc + [Fai Fc\FbFc) H^'lx- (A6) 

As will be shown below, when the gauge boson momenta are equal, which is the 
case in Eq. (|), the HTL-sums are related as 

= -^/fi^f- (A7) 

Inspection shows that the four-point function can be written in this case as 

^Tb = ^CUFa, Fb}^.:^'^,. (A8) 



The sum in Eq. ( |Al|) is most easily evaluated by applying the formula (see e.g. 



TT^nh) = dzf,{z)[F{z) + F{-z)], (A9) 

fco 

where rj = —1 corresponds to the Bose-Einstein and ?7 = +1 to the Fermi-Dirac 
distribution, the appropriate distribution chosen according to whether ko is even or 
odd multiple of irTi. The integration path encloses all the poles of the function F{z) 
on the right half of the z-plane. A calculation results in 



V / j = l j ijtj 



u;, = ^{k + p,y + ml (All) 

Next we perform a linear approximation in external momenta in the denominators and 
neglect the dependence of the external momenta elsewhere: 

1 " poo 



- n — - — -+n — \- (A12) 



Here K = (1, +k), K' = (1, — k), and the numerator /C is defined as ^ = K^^ ■ ■ ■ K'^"'^ 
(and similarly for /C'). Finally, after the /c-integral is performed we end up with 

m^->^"~^ = ^ / — ^l-vn ^ — - + -I]rT ^ — -1- (A13) 

6-2- J 4n \ j-J^) {P,- P^) . K 2j-,l\{P^-P^).K\ ^ ^ 

Here the sum over internal lines is divided to the fermion and boson contributions. 
Using the Eq. (|A13|) and plugging in the appropriate momentum configuration in 



the function gives 

n = fr + .?(T) / f . (A14) 

V J "^T^ K ■ PiK ■ P2J 

where Pi (P2)is the incoming (outgoing) momentum of the fermion and 

uJi{TY = C'p^ (A15) 



is the thermal mass of the left-handed fermion . 

The relevant momentum configuration of the function H^-i is Pi = P2 = 0, P3 = 
—Q and P4 = P, where Q is the momentum of the incoming and outgoing gauge boson 
and P is the momentum of the incoming and outgoing fermion. The term inside the 
braces in Eq. ( |A13|) may in this case be expressed as 

{...1 = --^ i . (A16) 

^ ^ 2{K-PyK-{P + Q) ^ ' 

On the other hand, the relevant momentum configuration of the function is Pi = 0, 
P2 = Q, P3 = P and P4 = P + Q, and the term inside the braces in Eq. (|A13|) becomes 



{...} = i-^ 1 + 1 . (A17) 

2{K-Py\K-iP-Q) K-{P + Q)j 

As the momentum Q is integrated over, and the propagator of the gauge boson is 
symmetric in Q, one can conclude from ( |A16| ) and ( |A17[ ) that H/i-i = — l/2if4;2. The 
four-point function has then the form 

rfr^ = u}{t){Fa,Fb] / . AI8 

>^ ^' at^{k-pyk-{p + q) ^ ' 

The effective propagators may be calculated along similar lines. The gauge boson 
vacuum polarization functions read as 

^^ABn^ =5^^3M(T)2(1-x2)(1-|l), 
Hf^ =5^^IIl = (5^^(^M(T)2 - Iht), (A19) 



where x = u/k and L = log((l + a;)/(l — x)). The thermal masses for the SU{2) and 
U{1) gauge bosons are up to the couphng constants equah 



m2 - Nf + m + 2 glT' 11 



Ml,,,, = ^^^^f^^ = ^^?r^ (A20) 

Here Nf is the number of left-handed SU (2) doublet fermion representations present in 
the thermal bath. We have performed the calculations with Nf = 12, i.e. the thermal 
bath is assumed to contain also top quarks. However, the results are not expected to 

9 



be very sensitive to small variations of Nf or 11./?,- The contribution in (|A20| ) 



come from the scalar doublet contribution. 

The propagator sum of the gauge bosons may be cast to the form 

^{F^, FsjAf J = + ^(r^'^r^'^ - Tl,)~\T^''' + r^'^tan^ 9^ + 2Tut^n9, 



(A21) 



where 



Ti 2 = uj -k - H/y - o 



12 



^?i^72(^) , (A22) 



and tan6'^ = gi/g2- 

The fermion propagators are of the form 

= i-P-^R-mM. (A23) 

where 

^lAP) = <r{T) l^§^ = '^ [Llo + (2 - xL)p . 7] , (A24) 

the notation being the same as in Eq. ( [A19| ), uji is the thermal mass from Eq. ( |A2(]| ) 
and a;^ defined in Eq. (^T[) gives the thermal mass of the right-handed neutrino arising 
from the scalar- fermion loop. 

Using Eqs. ( |AT^ ), ( |AT8D , QATOl) and (p^ , it is easy to deduce the following Ward 
identities: 



(Pi-P2Vr^(Pi,P2) = FA{MPir'-MP2)-'] 

g^g.r75(p,g) = {Fa, Fb} {mp + qy' - AiiP)-') . (A25) 

The latter of these identities holds only in the case where there is an integration over 
gauge boson momentum Q. 
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Figure captions 



Figure 1. Generic effective diagrams for the neutrino damping rates: a) diagrams 
contributing to the left-handed neutrinos, b) diagrams contributing to the right-handed 
neutrinos and c) non-contributing Diagrams. 

Figure 2. Damping rate of the left-handed neutrinos as a function of the vacuum 
expectation value of the doublet field: the total rate (solid line), contribution coming 
from particle- hole annihilation (dashed line) and the other contributions (dotted line) . 

Figure 3. Damping rates of right-handed hole-like excitations of neutrinos as function 
of the total scalar mass 777.5: (a) For Yukawa coupling Am = 0.1 and various ratios r = 
0, 2 and 6; (b) For fixed ratio r — 6 and for Yukawa couplings Am = 0.1, 0.5, 1. 

Figure 4. Damping rates of right-handed particle-like excitations of neutrinos as func- 
tion of the total scalar mass ms with ratios r = 0, 2: (a) Am = 0.1, (b) Am = 1-0. 



Table caption 

Table 1. The decay channels of the on-shell contribution to the damping rate of the 
right-handed particle (7+) and hole (7.) excitations and corresponding mass thresh- 
olds. 
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ms < uj'^ — uj^ 


ms > uj^ — uj^ 


ms > 2uj^ 


Up Svh 


VpVh S 
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ms> UJ^ - 


ms > 2uj~ 


UhS Up 


Vh^p S 





Fig. 1 a 





Fig. 1 b 



Fig. 1 c 
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